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Abstract 

Neutron Compton scattering measurements have the potential to pro- 
. vide direct information about atomic momentum distributions and adia- 

| batic energy surfaces in condensed matter. First applied to measuring the 

\Q . condensate fraction in superfluid helium, the technique has recently been 

extended to study a variety of classical and quantum liquids and solids. 
This article reviews the theoretical background for the interpretation of 
neutron Compton scattering, with emphasis on studies of solids. 

S. 

^ . 1 Introduction 

S3 
O 

Neutron Compton scattering is set apart from other branches of neutron scat- 
tering by the magnitude of energy and momentum transfers involved, typically 
in excess of 1 eV and 30 rA , respectively. Neutron scattering is renowned as 
a sensitive probe of collective properties in condensed matter, but in this ex- 
treme range of energies and momenta — the impulse limit — it is single-particle 
properties that are probed. The scattering occurs so rapidly, compared with 
the time scales of atomic motion in the sample, that the measured response is 
rather simply related to the equilibrium momentum distribution of the atoms. 

Thus, the term neutron Compton scattering refers, not to a distinct type 
of scattering, but to the usual neutron scattering cross-section in a limiting 
range of parameters. In the impulse limit, the scattered intensity, as a func- 
tion of energy for fixed momentum transfer Q, consists of a peak centred at 
E R = h 2 Q 2 /2M, the energy of recoil of a stationary nucleus on colliding with 
a neutron. The dependence of the peak position on the mass M of the struck 
atom implies that the scattering from different atoms appear at different ener- 
gies. This ability to separate contributions from different atomic species (and 
isotopes) is a useful feature of the technique. 

Doppler broadening results in a recoil peak whose width is proportional to 
the mean kinetic energy of the atoms, and whose detailed shape depends on the 
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distribution of atomic momenta. The principal tool in interpreting experiments 
is the impulse approximation, which predicts a precise relationship between the 
scattering data, in the form of the Compton profile, and the atomic momentum 
distribution. This gives an opportunity to compare results of experiments with 
predictions of realistic theoretical models of microscopic properties. 

As the name suggests, neutron Compton scattering is closely related to ex- 
perimental techniques in other branches of physics Compton scattering 
]2]||Q, i.e. inelastic scattering of X-rays, gamma rays or (10-60 keV) electrons 
provides information on electron momentum distributions, an example being 
the direct observation of the Fermi-Dirac distribution for conduction electrons 
in metals. In nuclear physics, an established technique for probing the momen- 
tum distribution of nucleons in nuclei is deep inelastic scattering |5||| of protons 
or electrons with energies of 100-1000 MeV and de Broglie wavelengths of the 
order of nuclear radii. At still higher energies (10 GeV and above), deep in- 
elastic scattering of electrons, muons or neutrinos fHHH probes the internal 
structure of nucleons and has been a key experiment in the confirmation of 
the existence of quarks. By analogy with these techniques, neutron Compton 
scattering is sometimes called deep inelastic neutron scattering. 

Neutron Compton scattering was first suggested |l(| as a technique for mea- 
suring of the condensate fraction in superfluid 4 He, and this has been an active 
field for three decades (see [jl.11.1^] and references therein). Recently, the de- 



velopment of spallation neutron sources, which have much higher intensities 
than reactor sources in the incident energy range required for neutron Comp- 
ton scattering, have opened up new applications. Recent experiments include 
studies of condensed noble gases [l^,14,l ^Jl^| , metals |17,[lq], normal liquid 4 He 
lHHliB lic l uid 3fIe 11 > solid 4 Hegf7 superfluid 3 He- 4 He mixtures @, 



and molecular hydrogen and deuterium |24,25,pqj2 



The extreme momentum transfers now experimentally reachable imply that, 
in some cases, the corrections to the impulse approximation are essentially neg- 
ligible. This suggests the possibility of direct and model-independent extraction 
of momentum distributions from Compton profiles, in contrast to the usual ap- 
proach of fitting theoretical predictions to the experimental profiles. Though its 
feasibility has not yet been demonstrated in practice, this approach, if success- 
ful, could establish neutron Compton scattering as the only technique capable of 
measuring atomic momentum distributions directly. A particularly interesting 
case is where the motion of the atom of interest is well described by an effective 
(Born-Oppenheimer) potential, or adiabatic energy surface. An example is a 
proton bound in a heavy lattice. In this case, the momentum distribution is 
the squared amplitude of the Fourier transform of the proton wave function, 



and from it, the potential energy function can be extracted [28,2 



The information obtainable in neutron Compton scattering is to some ex- 
tent complementary to that from diffraction experiments. The former measures 
the Fourier transform of a time averaged density; the latter the instantaneous 
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momentum density. Thus, for example, a neutron Compton scattering exper- 
iment on protons in a double-well potential of a hydrogen bond [2^,29] could, 



in principle, distinguish between a wave function with amplitude in both wells, 
and a statistical mixture of states in which the proton is localized in one well 
or the other. 

This article aims to summarize the theoretical framework for the interpreta- 
tion of neutron Compton scattering experiments, with an emphasis on applica- 
tions to solids. A thorough discussion of the physical principles underlying the 
impulse approximation is given in Sec. 2, based on the central concept of the 
scattering time. Sec. 3 addresses the problem of extracting information from 
the Compton profile. In the final section we discuss briefly the potential of the 
direct inversion approach. Although some examples are drawn from published 
experimental studies, no attempt is made here to review current experimental 
activity. For further background material, see Refs. pC|,31|. 



2 The impulse approximation 

We begin with an expression for the quantity which is measured in experi- 
ments, namely the partial differential cross section for the scattering of a beam 
of neutrons by a collection of atomic nuclei. Although the nucleon-nucleon 
interaction is very strong, it is sufficiently short-ranged that the scattering is 
a weak perturbation of the incident wave. Therefore, with the aid of a Fermi 
pseudo-potential, the scattering may be described in the first Born approxima- 
tion, and as a consequence the cross-section depends only on the changes in 
energy and wave vector of the neutron. We shall denote these by 



h 2 

huj = — {k 2 - k' 2 ) and Q = k - k', (1) 
2m 



respectively, where k and k' are the incident and scattered neutron wave vectors. 
A second simplifying assumption appropriate in the case of neutron Compton 
scattering is that the spatial scale of the scattering event, set by 1/Q, is too 
small to detect correlations between the positions of different nuclei, and hence 
that the scattering may be described to a good approximation as incoherent. 
This approximation is particularly good for scattering from protons, where 
the incoherent cross-section is larger than the coherent cross-section by almost 
two orders of magnitude. Under these circumstances the contribution to the 
scattering cross-section from nuclei of a particular species is proportional to the 
response function [32] 

1 r°° 

5i(Q,w) = 2^k y^^E^Q'*)' ( 2 ) 

where 

Yj(Q,t) = ( e -*Q-R, e *<W*)) (3) 
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is the density-density correlation function corresponding to the nucleus of atom 
j, whose position is represented by the quantum mechanical operator Rj. The 
corresponding Rj(t) has a time dependence in the Heisenberg representation, 

Rj.(t) = e im/h R je - iHt/h . (4) 

The angular brackets in Eq. (^) denote a thermal average of the enclosed expres- 
sion, as well as an implicit average over degrees of freedom which are passive in 
the scattering process, such as nuclear spin and neutron polarization states. 

The physical content of the response function can be made more transparent 
by passing to a somewhat less general expression. The thermal average in 
Eq. (|3|) involves, in principle, a sum over expectation values with respect to 
the complete set of stationary quantum states of the many-body Hamiltonian 
describing the nuclei and other particles, and their interactions. Let us suppose 
that the motion of a particular nucleus can be represented by an effective single- 
particle Hamiltonian which describes its interaction with its environment, and 
a corresponding set of single-particle states. We have in mind the case of a 
nucleus bound in a molecule, where an effective potential is constructed using 
the Born-Oppenheimer scheme. Let us denote the effective single-particle states 
by |n) and their energies by E n . The contribution to the response function from 
a single nucleus reduces to 

Si(Q,u) = ]T Z- V^/^IHe-^V)! 2 ^ " (En' ~ E n )), (5) 

nn' 

where Z is a normalization factor (the thermodynamic partition function). In 
this expression we may recognize a sum over transitions from initial states re, 
weighted by a thermal Boltzmann factor, to final states n', of a transition prob- 
ability. The latter is given by Fermi's golden rule as the product of a squared 
matrix element and an energy-conserving delta function. The scattering is rep- 
resented by the operator e _i< ^ R which couples the plane wave of the neutron 
with the position of the nucleus. 

Returning now to the general expression, Eq. (||), the form of the response 
function as a sum of separate contributions from each nucleus reflects the nature 
of the incoherent approximation, which neglects correlations between different 
nuclei. The impulse approximation, which is the main topic of this section, 
consists in a similar neglect of time correlations in the motion. To be specific, 
if the time scale of the scattering event is much shorter than the characteristic 
time of atomic motions in the sample, then the nuclei may be regarded as free 
particles in so far as the scattering probability is concerned. As a result, the 
latter depends only on the momentum of the nucleus in its initial state. Crudely 
speaking, then, we may summarize the approximations appropriate for neutron 
Compton scattering as follows: the Born approximation, that each neutron 
scatters only once; the incoherent approximation, that each scattering process 
involves only one nucleus; and the impulse approximation, that the neutron 
remains in the vicinity of the nucleus for a time too short to sense anything 
except how fast the nucleus is moving. This last approximation is the most 
subtle of the three, and it will be discussed in detail now. 
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2.1 The short time expansion 



Motivated by the idea that only short times are relevant, let us suppose that 
the correlation function Y(Q,t) for a single nucleus is dominated by its be- 
haviour for small t. The time dependent Heisenberg operator R(t) has a Taylor 
expansion, 

R(t) = R + (P/M )t + y 2 (F/M)t 2 + ..., (6) 

where P = MR = iM[H, R]/7i is the momentum of the struck nucleus which 
has mass M, and F is the force defined similarly. All the operators in the 
expansion are evaluated at t = 0. We remark that the identification of MR 
with momentum is valid only in the absence of velocity-dependent forces, i.e. it 
is assumed that there are no magnetic fields present, and that the motion of 
the nucleus is non-relativistic. 

Reserving for later the question of the precise criteria for validity of the 
present mathematical procedure, let us examine the form of the response func- 
tion resulting from the neglect of terms of order t 2 and higher. This amounts 
to an assumption that the struck particle is effectively free so that F = 0. With 
the aid of the operator identity 

e A+B = e A e B e -y 2 [A,B ] ^ (7) 

which holds when [A, B] commutes with both A and B, the correlation function 
is found to reduce to 

Y(Q,t) = ( e i(Q-P/A/+^fl)t) ) ( 8 ) 

where lor = HQ 2 /2M is the (free atom) recoil frequency. The exponent has an 
interpretation in terms of the kinematics of a neutron-nucleus collision. If the 
momentum of the nucleus is P before the collision, then it is P + HQ after the 
collision, and energy conservation requires 

uj = [(P + HQ) 2 - P 2 ]/2Mh = Q • P/M + lor. (9) 

In particular, hujR is the energy imparted to a stationary nucleus by a collision 
with a neutron. 

If the nucleus were indeed at rest before the scattering, Eq. (|8|) would give 
Y(Q,i) = e lujRt , and the response function would consist of a delta function 
line at the recoil frequency. In fact, the nucleus will be in a quantum state 
having a distribution of initial momenta, and the line will therefore be Doppler 
broadened. Each possible initial momentum results in a contribution to the 
scattering intensity at a frequency given by the conservation constraint, Eq. (||). 
Specifically, defining the momentum distribution by 

7»(q) = <$(q-P/ft)>, (10) 
the response function in the impulse approximation is 

5i(Q, to) = h- 1 [ n(q)£(w - hQ ■ q/M - u R ) dq. (11) 
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This result, that the scattering cross-section is directly related to the momentum 
distribution of the struck nuclei, is of central importance in neutron Compton 
scattering experiments. 

It is important to emphasize that, although it represents a single-particle 
response, the momentum distribution n(q) is actually a property of the many- 
body system of all the nuclei and their interactions. In other words, the momen- 
tum distribution of a single nucleus depends on its environment and therefore 
on the behaviour of the system as a whole. Let us examine again the special 
case where the motion of a nucleus is given by a set of effective single-particle 
states |n). Here, Eq. (|i~0|) becomes 



(12) 



n(q) = {l/2^fY, z ~ le ~ En/kBT J e^ q r (r|n)dr 



At zero temperature, this reduces to the square modulus of the Fourier trans- 
form of the ground state wave function. Thus, there is a direct relationship 
between the measured quantity, Si(Q,uj), and the quantum mechanical wave 
functions of the nuclei in the sample. 



2.2 The scattering time 

Looking at the short time expansion of the nucleus' position operator, Eq. (||), 
we see that the first term neglected in the derivation of the impulse approxi- 
mation is proportional to the operator F representing the force experienced by 
the nucleus. Denoting root mean square values by an overbar, the corrections 
are expected to be small if 

F^Qr s «FQ, (13) 

where r s is a quantity with the dimensions of time, identified as the time scale 
of the scattering process, or scattering time. At first sight, it might appear 
reasonable to relate t s to the time taken for a neutron wave packet to pass 
the vicinity of the nucleus. The corollary, that the degree of coherence of the 
neutron beam plays a role in deciding the validity of the impulse approximation, 
is largely erroneous, as will be demonstrated presently. 

To establish the significance of the scattering time, let us adopt the op- 
erational point of view that, since the aim is to approximate the correlation 
function Y(Q, t), the relevant time scale should be obtained from Y(Q, t) itself. 
The response function Si(Q,uj), as we have seen, consists in general of a peak 
centred at the recoil frequency ojr, with a certain width Aw. It follows that 
the general form of Y(Q,t) is of an oscillatory function e lu)Rt , modulated by a 
decreasing envelope of width 1/Au (Fig. 1). The fact that Y(Q,t) goes to zero 
confirms our expectation from physical arguments that knowledge of the struck 
atom's motion for a limited time span is sufficient to predict the scattering 
response. Accordingly, we identify the scattering time t s as the time scale for 
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Figure 1: Schematic form of the correlation function Y(Q, t), whose Fourier transform 
is proportional to the neutron scattering cross section. The frequency of the oscillations 
is the reciprocal of the recoil frequency of the scatterer, and the amplitude falls off on 
a scale of the scattering time r s . 



the decay of the correlation function Y(Q,i) to zero, which equals 1/Au, the 
reciprocal of the recoil peak width |jl2|| . If it happens that Si(Q,uj) is highly 
structured, i.e. has features on several frequency scales, this means that Y(Q, t) 
depends on more than one characteristic time. The overall decay envelope is 
determined by the longest of these characteristic times, and it follows that t s 
is the reciprocal of the width of the narrowest feature in Si(Q,u). 

Strictly speaking, for scattering from solids the correlation function Y(Q,t) 
does not go to zero, but to |(e l Q' R }| 2 , corresponding to the amplitude of the 
elastic line in Si(Q,u) |32|. However, this constant component does not affect 
the argument: since we are interested in the scattering at energies near oj = ujr, 
we may discard the elastic scattering and take Y(Q,t) to consist only of the 
inelastic part, which tends to zero as described. 

With t s defined in terms of the structure of the response function, a rigorous 
determination requires, in principle, a detailed calculation of Si(Q,co) for the 
system of interest. In the absence of a detailed theory, some rough estimates 
can be made. A convenient procedure is to estimate r s self-consistently using 
the impulse approximation itself. We take for Si(Q,u) a peak centred on lor 



of width Aw. According to Eq. (11), Aoj is proportional to the width of the 
distribution of Pq, the projection of momentum P along the direction of Q, 
and thus 

M . . 

T ' ~ WW' ( ' 

The inverse dependence of the scattering time on Q is an essential feature. It 
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implies that, as anticipated on physical grounds, the impulse approximation 
becomes exact in the limit of large wave vector transfer. This statement is true 
provided the forces on the nuclei are always finite, for if F can grow arbitrarily 
large the condition in Eq. (|l^) is not satisfied no matter how short the scattering 
time. We shall therefore not consider pathological cases, such as scattering 
from a "hard-core" fluid (3^] or from a particle undergoing Brownian motion 
]34|| , where corrections to the impulse approximation do not become negligible 
as Q — > oo. 



With our estimate, Eq. (14), for the scattering time, the criterion for validity 



of the impulse approximation reads 

F^d/Q) « (P%)/M, (15) 

where Fq = F • Q/Q is the projection of the force onto the direction of Q. This 
has the following interpretation. If 1/Q is regarded as the length scale of the 
scattering event, it is required that the work done by the force Fq in moving 
the nucleus this distance be negligible compared with the root mean square 
kinetic energy of the nucleus due to motion along Q; in other words, forces on 
the nucleus do not change its energy appreciably during the scattering process. 
This feature of impulsive scattering is reflected in the energy-conserving delta 



function in Eq. (11), which would otherwise involve a contribution from the 



nucleus' change in potential energy. 

An alternative discussion of the short time expansion can be made within 
the framework of the Gaussian approximation for incoherent scattering [p5| , 
in which the correlation function is written Y(Q,t) = exp [-Q 2 w(t)\, with the 
width function w(t) proportional to the mean square displacement of the atom 
at time t. The impulse approximation is obtained in an asymptotic analysis by 
expansion of w(t) about a saddle point in the complex plane. The correction 
terms may be estimated in terms of the moments of Si(Q,u), and are small 
if the "skewness" of the recoil peak, related to the third moment, is small. 
The condition for validity of the impulse approximation obtained in this way is 
equivalent to Eq. (|l5|). 



2.3 Bound nuclei 



The preceding derivation of the conditions for validity of the impulse approxi- 
mation is not rigorous, in that the use of the impulse approximation to predict 
its own range of validity is a circular argument. An example in which the rea- 
soning fails is that of scattering from a nucleus bound in a harmonic potential. 
This example will now be treated in some detail. It will lead us to extend 
somewhat the concept of scattering time, and will serve as an introduction to 
a general discussion of scattering from nuclei bound in solids. 

For the harmonic oscillator, 5j(Q,u;) can be calculated without approxima- 



tion [32]. It will be sufficient for our purposes to examine the case of zero 
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Figure 2: Solid vertical lines: the response function for scattering from a nucleus 
bound in an ideal harmonic well of frequency ujq (each delta function contribution is 
plotted as a vertical line of height equal to the corresponding weight). The momentum 
transfer Q is such that the recoil frequency con is 15luq. Dashed curve: the same 
quantity as predicted by the impulse approximation. 



temperature, for which the exact result is 

oo 

Si(Q, u) = h^e-"*/" £ [(u R /u; ) n /n\\ 5(u - tiuj q ). (16) 

n=0 

The prediction of the impulse approximation is 

Sj A = (2irh 2 oj oJ R )~ 1/2 exp 

and the self-consistent estimate of its range of validity, from Eq. (|l5|), is u R ^> 
ujq. The two functions are compared in Fig. 2 for uj r /ojq = 15. The exact 
Si(Q, to) is an array of infinitely sharp lines, whereas the impulse approximation 
result is a smooth Gaussian, drawn as a dashed curve. The impulse approxi- 
mation fails spectacularly to reproduce the fine structure of the spectrum. 

As is evident from the general expression, Eq. (||), for the response function 
in terms of single particle states, the source of the problem is energy quantiza- 
tion. Each sharp line represents the absorption of a number of quanta by the 
oscillator. Since Si(Q,uj) contains infinitely narrow structure, the scattering 
time, as defined above, is infinite, and the short time expansion fails. 

Not all is lost, however. Evidently, the impulse approximation does give 
an accurate account of the envelope of the palisade of delta functions. The 



(oj - io R f 
2u Q u) R 



(17) 
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agreement improves on increasing Q. In other words, the impulse approximation 
describes 5j(Q,u;) on a frequency scale large compared with uq. In practice, 
such a coarse description is likely to be adequate. For nuclei bound in molecules 
or solids, even if broadening mechanisms intrinsic to the system (discussed 
below) are insufficient to smear out structure on a scale of ujq, achieving the 
instrumental resolution needed to discern the separate lines would most likely 
be a difficult task, with an incident energy large enough to accomplish the 
experiment. 

The ability of the impulse approximation to describe the envelope of the 
discrete spectrum is not restricted to the ideal oscillator. The generalization to 
arbitrary potential wells is straightforward j3(|. At zero temperature, Eq. (||) 
is 

Si(Q,uj) = ]T |<0|e- iQ - R |n)| 2 ^ - (E n - E )). (18) 

n 

The scattering intensity is concentrated at values of u = {E n — Eq)/T% for which 
the matrix element is large. Now, the essence of the impulse approximation is 
that the momentum transfer, and hence the energy transfer, is so large that 
the nucleus' final state is approximately that of a free particle, a plane wave. 



Substituting a final state e iq R with E n = h 2 q' 2 /2M into Eq. @ yields |7| 



Si(Q,w) = J n(q)5[hiJ-h 2 (Q + c l } 2 /2M + E ]dq, (19) 



where re(q) is the momentum distribution defined previously in Eq. fll2j). For 
each q there is a contribution to Si(Q,uj) at a position relative to the recoil 
frequency 

oj -UR = hQ-ci/M + hq 2 /2M - E Q /H. (20) 

For sufficiently large Q, the first term dominates the others on the right hand 
side, and hence Eq. (|l^) is essentially identical to the impulse approximation 
derived using the short time expansion. 

In this re-derivation of the impulse approximation, three approximations 
were made. The first was the replacement of the final state by a plane wave in 
the evaluation of the matrix element in Eq. (18). Clearly it is sufficient that 



this approximation be accurate in the region of space where the ground state 
wave function is appreciable, near the centre of the potential well. That this 



is indeed the case follows from the ideas of the WKB approximation [38|, in 
which the wave function of a bound state is written as a standing wave with a 
local wave vector q'(x) = [2M(E n — V(x))/h 2 ] 1 ^ 2 . The impulse limit requires 

V(x) « E n (21) 

within the spatial extent of the ground state, since then q'{x) is approximately 
independent of position. In addition, to ensure that the WKB wave function is 
a reasonable approximation of the excited state we require the fractional change 



in q'{x) over one wavelength 27r/q' to be small [38|. In the present context this 



requires dV/dx <C q'E n , which is a weaker constraint than Eq. (p|) 
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The second approximation made in re-deriving the impulse approximation 
was ignoring the discreteness of the spectrum of final states. That is, the sum in 
Eq. ( Jl8| ) was replaced by an integral over all momenta. It is here that the delta 
function structure of Sj(Q, a;) is lost, and the procedure is therefore justified if 
our aim is to calculate the envelope of the spectrum. 

The final approximation was the neglect of the second and third terms of 
Eq. (20). Since E -h 2 q 2 /2M is, on average, the potential energy in the ground 



state, we obtain the condition 

(V) « hQ{P 2 )^ 2 /M. (22) 



This is evidently a stronger condition than Eq. fl21|) since the right hand side 
increases only linearly with Q. Hence, if Eq. ( |2^ ) is satisfied, all the approxi- 
mations leading to the impulse approximation are justified. 

The condition may be expressed in terms of the forces on the atom by not- 
ing that F = — VV. If F represents the root mean square force in the ground 
state, which has a spatial extent of order %/ (P 2 ) 1 / 2 , then the ground state po- 
tential energy is approximately HF / (P 2 ) 1 ^ 2 . Inserting this estimate in Eq. (^) 
and ignoring any dependence on the direction of Q yields F <C Q(P 2 )/M, in 



agreement with the self-consistent estimate, Eq. (15) 



We conclude that the short time expansion, and the associated self-consistent 
assessment of its validity, are adequate for scattering from a bound nucleus, 
provided we are satisfied with a description of the envelope of Si(Q, u). In this 
context, it is appropriate to define the scattering time r s as the inverse width 
of the narrowest envelope, which is the time scale of processes during scattering 
which determine the broad structure of the response function. 

The present discussion has considered atoms bound in potentials with in- 
finitely high walls. In practice, of course, a sufficiently large impulse will eject 
the atom from the molecule or lattice which binds it. It has been suggested 
in the literature that the impulse approximation is valid only when the recoil 
energy is large compared with the binding energy Eb- This is indeed the con- 
dition necessary for the scattering response to be a smooth function as the im- 
pulse approximation predicts, rather than an array of narrow lines. However, 
the conclusion of this section is that there is an intermediate energy range, 
ujq <C ujfi < Eg/h, in which the impulse approximation provides an accurate 
description of the envelope of the response function. 



24 Solids 

We turn now to the subject of scattering from nuclei bound in a lattice, in which 
motions of different nuclei are coupled. In other words, we wish to generalize the 
preceding discussion, which treated lattice vibrations as a collection of indepen- 
dent oscillators (i.e. in an Einstein model), to a situation where the vibrational 
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modes are collective in nature. We aim to provide qualitative estimates for 
the range of Q for which the impulse approximation is an accurate description 
of the response function. The task of applying the impulse approximation in 
detailed calculations for particular model systems will not be attempted here; 
see, for example Refs. [S3] and |39[|. 



The first point to be made is that the picture of Sj(Q, u) as an array of delta 
function lines is no longer appropriate when atomic vibrations are coupled. The 
reason is simply that final states excited by the scattering are drawn from a 
continuous distribution, namely the phonon density of states, rather than from 
a set of quantized levels as for the oscillator. To be specific, let us consider 



a perfect harmonic lattice, at zero temperature. Referring to Eq. (18), the 
final states are now taken to be excited states of the lattice as a whole, rather 
than single-particle states. It is convenient to group together terms with a 
given number of phonons excited. The leading term, then, has u = and 
yields a delta function elastic line just as for the harmonic oscillator. The one- 
phonon contribution, however, is not a sharp line but a continuous distribution 
proportional to the phonon density of states Z(u). Explicitly p2| , |3S| ] , 

$(Q,a;) = e~ 2W ^ {5(lo) + (u R /u)Z(u) + ...}, (23) 

where e~ 2W ^ is the Debye-Waller factor. Continuing the series, the two- 
phonon term is the product of a matrix element e~ 2W {ujr/uj) 2 and the two- 
phonon density of states, 

POO 

/ Z(lo')Z(lo - to') dw', (24) 
Jo 

which is the convolution of Z with itself. Higher terms in the phonon expansion 
are proportional to repeated self-convolutions of Z, and therefore become pro- 
gressively broader and smoother. If Q is large, by the time energy transfers of 
the order of the recoil energy are reached, the number of phonons excited is so 
large that the central limit theorem applies and the n-phonon density of states 
is essentially Gaussian. The different contributions merge to form a smooth 
recoil peak. These features are illustrated in calculations by Mayers, Andreani 



and Baciocco [ 40 1 using a Debye model, and by Evans et al. [41 and Fielding 
et al. p|] using realistic densities of states for ZrH2 and ZrD2- 



The smoothing effect of the continuous distribution of phonon states is a "de- 
phasing" effect, i.e. Y(Q,t) decays because the spread in frequencies causes a 
loss of coherence over time. This is to be contrasted with decay of correlations 
due to damping mechanisms. The latter are a consequence of interactions of 
phonons with impurities, with electrons, and with each other (due to anhar- 
monicity), resulting in a finite phonon lifetime. Such effects may be included in 



Eq. (18) by replacing the delta function by a peak of width equal to the inverse 



phonon lifetime. 

The effects of finite temperature are to broaden the scattering response still 
further. They enter through both de-phasing and damping mechanisms: the 
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former because of the ksT spread in initial energies, the latter because the 
phonon lifetime is a decreasing function of temperature. 



Let us summarize the various energy scales that have arisen in our discussion 
so far. The characteristic frequency of the phonons is the Debye frequency ljd, 
corresponding to an energy scale typically of the order of a few hundred kelvin. 
For low temperatures, ksE <S Tiujdi the mean kinetic energy per atom, Ek, is 
of order 3?iW£)/4, while for high temperatures it approaches the classical value 
3/cbT/2. The lifetime r of phonons is related to the thermal conductivity of the 
solid by K = Cv 2 /3t, where C is the heat capacity and v is the sound veloc- 
ity EJ. Inserting representative values leads to an estimate that the lifetime 
broadening fijr is typically 1 K or less, and hence negligible in comparison to 
other energy scales in the system. 

The energy scales related to the scattering process are the recoil energy Twjr, 
typically in excess of 1 eV, and the width Auj of the recoil peak. These are the 
scales which determine the applicability of the impulse approximation. In view 
of the smoothing effect of the continuous phonon density of states, we judge 
that the self-consistent criterion for the validity of the impulse approximation 
is appropriate. In order to apply it, we require an estimate of the forces on an 
atom. For a harmonic oscillator, the mean square of one Cartesian component 
of the force is M 2 ojq{x 2 ) = 2Mu>qEk /3, and thus for a harmonic or nearly 
harmonic solid Fq is of the order of (MajfjEx) 1 ^ 2 ■ The criterion for validity of 
the impulse approximation becomes 

Q » (Mu 2 D /E K ) 1 ' 2 ~ {Muj D /h) 1 ' 2 , (25) 

where the second estimate is valid for temperatures low compared with the 
Debye temperature. The quantity on the right is the inverse of the root mean 
square displacement of a harmonic oscillator of frequency cj_d/2. Thus, the 
criterion for the impulse approximation is more easily satisfied if the atoms are 
weakly bound, as one would expect. As remarked previously [|^J|^], Debye 
frequencies do not depend strongly on atomic mass, and hence the Q value 
required to reach the impulse limit should increase with mass as M 1 / 2 . 

Neutron Compton scattering investigations of monatomic solids reported 



recently include experiments on argon, krypton and xenon [13, lithium [17|, 



sodium [ 18 1 and 4 He [21]. To take a representative example, lithium [17| has 
a mass of 7 a.m.u. and a Debye temperature of approximately 400 K, yielding 
the criterion Q S> 8 rA _1 for the impulse limit. For Q ~ 100 rA _1 , the highest 
used in the experiment cited, the recoil frequency ujh ~ 3 eV, and the recoil 
peak width Aui is about 0.3 eV. 

These considerations are readily generalized to more complicated systems. 
An example is the hydrogen molecule, which has been the subject of recent 
experiments [24,^5,2^,27]. In this system the binding of protons within an H2 



molecule is much stronger than the forces between molecules. The mismatch 
of energy scales is reflected in the spectrum of vibrational states, which is en- 
visaged as consisting of two bands: a broad acoustic band, corresponding to 
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molecular motions characterized by the Debye frequency up, and a narrow high 
frequency band centred at the intramolecular vibrational frequency uy. The 
width of the narrow band is of the order of Wp/wy. In H2, oJy is nearly two 
orders of magnitude greater than cjd, and the molecular and intramolecular 
motions are effectively decoupled. The characteristics of the scattering are thus 



dominated by the intramolecular vibrations, and in fact Mayers [24] found the 
scattering from liquid H2 to be indistinguishable from that from the polycrys- 
talline solid. In the low temperature limit (ksT <C hujy) appropriate here, the 
criterion for validity of the impulse approximation is Q ^> (Muy /h) 1 ^ 2 ~ 10 
rA _1 , well within the experimental range. 

Similar arguments apply for scattering from light atoms in a heavy lattice 
|44], such as the proton in KHCO3 [45|. Here the proton vibrational modes are 



found in a high frequency band, narrower than the acoustic band by approxi- 
mately (M1/M2) 1 / 2 , and are effectively decoupled from other motions because 
of the large mass ratio. 



2.5 Liquids 



The case of scattering from liquids is the most well developed and, at the same 
time, the most controversial, application of the impulse approximation. Here 
we restrict ourselves to a few general remarks; details may be found elsewhere 

I Tgjll|,0,|8|,||j5^f3ilil'lll'lll'lll'il 1 • 



For a monatomic classical liquid, such as a condensed heavy noble gas |13[ , 
the application of the impulse approximation appears straightforward. For ex- 
ample, taking order of magnitude values for the force between atoms in the 
liquid near its triple point, (F 2 ) 1 / 2 ~ 10 3 KrA -1 |l3|], and the kinetic energy 
Ek = 3/c#T/2 in the range 200 to 400 K, we obtain an estimate Q > 5 rA _1 
for the impulse limit to be reached. Indeed, in the experiments cited, devia- 
tions from the impulse approximation were observed to be small at momentum 
transfers between 17 and 29 rA _1 . 



The situation appears to be not very different for most monatomic quantum 
fluids, except, of course, that the evaluation of the kinetic energy must take 
account of the quantum zero point motion of the atoms. Examples of recent 
experiments on normal liquids in which quantum effects are important include 
studies of 4 He [piraj2qM| , 3 He pj and neon MM. 



One way to estimate orders of magnitude is to use a "cell model" of the 
liquid, in which, on short time scales, an atom is assumed to move in a roughly 
spherical cage created by neighbouring atoms, with intermolecular potentials 
of, say, Lennard-Jones form. For example, Andreani et al. ]^0| have argued that 
such a model, with the total interatomic potential represented by an effective 
harmonic vibrational frequency ujq = 14 K, accounts reasonably well for the 
observed temperature dependence of the mean kinetic energy of normal liquid 
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4 He. 

Assuming that the same model suffices to estimate the forces on an atom, 
we find that Q S> {Muq/Ti) 1 / 2 ~ 1 rA -1 is the condition for the impulse ap- 
proximation to be accurate. In fact, this type of estimate is highly misleading. 
The impression that the impulse limit is attained for rather low Q results from 
the low mass and the fact that the attractive part of the Lennard-Jones poten- 
tial is very weakly binding. However, the experience gained from a variety of 



theoretical and experimental investigations over the past few decades (see [11 
and references therein) has shown that the repulsive "hard-core" -like part of 
the interatomic potential is crucial in determining the validity of the impulse 
approximation. Indeed, as mentioned previously, for a true hard-core fluid the 
impulse limit is not reached no matter how large the momentum transfer |3 



For the most extreme example of a quantum liquid, 4 He in its normal and 



superfluid phases, the work of Silver |p4,5q,5q| is a definitive theoretical study of 



corrections to the impulse approximation. Its conclusion is that, although the 
correction terms in a formal expansion of Si(Q,uj) as Q — > oo are proportional 
to powers of 1/Q |^9| , the nearly hard-core nature of the interactions results in a 
broad range of crossover to the asymptotic limit, in which the corrections scale 
as log Q. As a result, the corrections are not negligible, even for the highest 
Q values that might conceivably be attained in experiments. Silver's work also 
provides a systematic method of calculating the corrections, which has been 
applied successfully in measurements of the condensate fraction in superfluid 
helium [57]. 



For non-monatomic liquids, there are complicating (and interesting) features 
arising from the internal structure of the molecule. Recent work includes ex- 



periments on molecular hydrogen and deuterium [ 24, 25, 26, 27 1. 



2.6 Corrections to the impulse approximation 

Deviations of the response function Si(Q,uj) from the prediction of the impulse 
approximation are often collectively termed "final state effects", since they re- 
sult from deviations of the final state of the scattering process from plane wave 
form. In fact, since Si(Q,u) reflects properties of initial as well as final states, 



and deviations arise from both sources [37.4C], we shall refer simply to "correc- 
tions to the impulse approximation" . 

In general, the corrections take the form of a broadening of the response 
function; this is in part an effect of a finite lifetime of the final state, due to 



collisions of the struck atom with its neighbours |1C[. An example is the broad- 
ening effect of the phonon lifetime discussed in Sec. 2.4. The corrections are 
frequently embodied in a "final state resolution function", which is convolved 
with the result of the impulse approximation to obtain the predicted Compton 
profile. 
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When deviations from the impulse approximation are appreciable, accurate 
calculations of the final state resolution function require a detailed theory of 
interactions in the system under study. An example is the hard-core pertur- 
bation theory [54|,55,5(| for normal and superfluid 4 He, which is a quantum 



many-body theory of the fluid incorporating scattering data for the helium in- 
teratomic potential. Here, we shall not enter into details of such calculations, 
but make a few general comments, principally aimed at the case of scattering 
from solids for which the corrections are small. 



Let us recall the approach taken in Sec. 2J5 and Ref. 36], where the impulse 
approximation was derived under three assumptions: that the final state is 
nearly a plane wave, that the discreteness of the final state energies is unimpor- 
tant, and that the last two terms in the energy conservation condition, Eq. (pC|), 
may be neglected. The last assumption, that 

E m h 2 q 2 /2M, (26) 



is exactly true only for a free particle system. Mayers [37] has termed the de- 



viations from Eq. (26) "initial state effects". Here we focus on the fact that, 
as mentioned previously, Eq. (^6|) is the strongest of the three approximations 
made, and therefore if we relax it, the result is likely to be a better approxima- 
tion than the impulse approximation. For example, if the final state energy E n 
is set equal to its WKB estimate % 2 q 2 jIM + V(x), and if Eq — V(x) is replaced 
with its average value, namely the mean kinetic energy Ex, one obtains 

Si(Q, w) = J n(q)5[huj - h 2 (Q + q) 2 /2M + E K ) dq, (27) 



a result first suggested by Stringari |5q[ . It is a better approximation than the 
impulse approximation at low temperatures, but becomes less useful at higher 
temperature, since the replacement of a distribution of energies by an average 
ceases to be valid EL 



The Stringari formula shows that one aspect of the corrections to the impulse 
approximation is a shift of the maximum of the recoil peak to lower frequency. 
Such a shift is indeed visible in Fig. 2, for example. In addition, the correc- 
tions induce an asymmetry in the peak. This suggests that a straightforward 
symmetrization of the data can be used to remove partially deviations from 
the impulse approximation. In the Sears [49] method, symmetrization is the 



first step in a systematic self-consistent correction procedure. The method is 
based on a formal asymptotic expansion of S*j(Q, uj) in powers of 1/Q, in terms 
of successive moments of the final state resolution function. It is found that 
symmetrization of the measured recoil peak eliminates corrections of order Q _1 , 
leaving residual terms of order Q~ 2 . The antisymmetric part of Si(Q,uj) is pro- 
portional to Q^ 1 and amounts to a measurement of deviations from the impulse 
approximation; this information can then be used to correct for the residual de- 
viations in the symmetrized data, and the result is a response function corrected 
up to order Q~ 2 . Comparisons of correction methods for experimental data on 
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bound light atoms have been made by Evans et al. Kl[ and by Fielding et 
al. @. 



An analysis method recently developed by Glyde [^] , also based on moment 
expansions, aims to extract both the limiting impulse approximation response 
function, and measurements of the finite Q corrections, from experimental neu- 
tron Compton scattering data. In this approach, rather than concentrating on 
extreme momentum transfers in order to minimize deviations from the impulse 
approximation, one collects data over a broad range of Q. It is the measurable 
difference in the Q dependence of various contributions to the expansion mo- 
ments that allows the Q — > oo (impulse approximation) component of the data 
to be isolated, and the corrections to be estimated. This technique has been 
applied in a recent study of normal liquid 4 He and liquid neon p6| . 



3 The Compton profile 

In this section we turn to the problem of extracting information about the 



momentum distribution from measured data. We shall take the approach [28,B£] 
of assuming that the impulse approximation is accurate. As we have seen, 
this assumption is justified in the case of scattering from solids, where the 
impulse limit is well within the experimental range of momentum transfers and 
where leading order "final state" corrections can be handled, for example, by 
symmetrization of the data. 

The measured response function Si(Q,u>) is given in the impulse approxi- 
mation by Eq. (|TT|), Let us define J = (fr 2 Q/M)Si and rearrange the delta 
function to give 

J = J n(q)<f(q • Q - y) dq, (28) 

where Q is the unit vector along Q, and 

y={M/%Q){uj-uj R ). (29) 

This form of Si(Q,Lu) has an important consequence. Consider an isotropic 
system, where the scattering is independent of the direction of Q. Then J 
depends only on the variable y, rather than on Q and uj independently. In 
other words, the data "collapse" onto a function J(y) of a single variable. This 
phenomenon is known as y-scaling, and J(y) is termed the Compton profile. It 
is seen from the delta function in Eq. (|2g| ) that y is the projection of the atom's 
momentum onto the scattering vector. 

The function J(y) is a convenient and standard form for experimental results 
to be presented. The recoil peak is shifted to be centred at y = 0, and if the 
profile is normalized to unity the mean kinetic energy per atom is directly 
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related to the second moment of the Compton profile by 



/oo 
y 2 J{y)dy. (30) 
-oo 

This result follows from the general moment relations for the incoherent re- 
sponse function [32]. 

For an anisotropic system such as a single crystal, the response function 
depends on the direction of Q as well as on y, and one then speaks of the 
directional Compton profile, J(Q,y). It has a second moment related to the 
kinetic energy associated with motion along Q. 



3.1 Reconstructing momentum distributions 

The Compton profile and the momentum distribution are related by Eq. 
Consider J(Q, y) as a function of y for fixed Q, and define a coordinate system 
xyz such that Q points along the z-axis. Then 

J{Q,y) = J n(q x ,q y ,y)dq x dq y . (31) 

The Compton profile is seen here to be the integral of n(q) over a plane normal 
to the vector Q, i.e. a projection of the momentum density onto the direction 
of the momentum transfer. This mathematical relationship between J(Q,y) 



and n(q) is known as a Radon transform [6C|. An important property of the 
transform is that it is invertible: given the directional Compton profile for all 
values of its arguments, the momentum distribution can be extracted. 

The Radon inversion formula can be expressed in many mathematically 
equivalent forms, but these are not equivalent in practice, when the data are 
finite and affected by noise and instrumental resolution. The reconstruction 
problem for the Radon transform, and the associated mathematical questions 
of stability, uniqueness, accuracy and resolution, have been studied thoroughly 
in connection with computerized tomography in diagnostic radiology pl| . We 
note, in passing, that the Radon transform can be expressed in the framework 
of wavelet theory |p2[ , but the consequences, if any, for the reconstruction prob- 
lem do not appear to have been explored. Here we describe a reconstruction 
method proposed by Reiter [28,29| based on the work of Davison and Grun- 
baum [63 134 1, which is a variant of the "filtered back-projection" technique in 
common use in medical applications of tomography. 

The method involves decomposing the angular (Q) dependence of the data 
into spherical harmonics, and the y-dependence into products of Gaussian and 
Hermite functions. Specifically, 

J(Q, y) = vr-VSg-^,,) 2 £ A nlm H 2n+l (y/y )Y lm (Q), (32) 

nlm 
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where Yi m are the usual spherical harmonics, as defined, for example, in Ref. 
The constant yo is chosen to match the width of the Compton profile. As may 
be proved by direct integration [29[|, this series is the Radon transform of 



n(q) = vr 3 /V^) 2 ]T B nlm {q/y ) 1 L l + l /\q 2 /j/g)n»(q), (33) 



nlm 

where 



Bnlm = (-l) n 2 In+i n\A nlm . (34) 

Here H and L are Hermite and Laguerre polynomials, with the standard nor- 
malization |35| that the coefficient of x m in H m (x) is 2 n , and the coefficient of 
x m in L c f n {x) is (— l) n /nl. The method is thus to determine the {^4} coefficients 
from the data, to generate the {B} coefficients using Eq. (|34"|), and finally to 
obtain the reconstructed momentum distribution in the series form in Eq. ([^) . 

Two properties of this inversion method make it particularly suitable for 
application in neutron Compton scattering. The first is that instrumental reso- 
lution is readily accounted for by including the resolution function in the fitting 
process used to extract the A coefficients. In other words, the functions actu- 
ally used to fit the data are the convolutions of those on the right hand side 
of Eq. p^ ) with the instrumental broadening function. The second desirable 
property is that the expansion functions are very well matched to the charac- 
teristics of the data. Indeed, if the underlying binding potential is isotropic 
and harmonic, then J(Q, y) is an isotropic Gaussian and only the first term, 
proportional to Aqqq, is needed. Even for anharmonic potentials, the Compton 
profile is expected to consist of a compact peak, and a fairly small number of 
terms in the expansion should suffice. The {^4} coefficients, therefore, are an 
economical description of the data, and the coefficients with n > are directly 



related to the anharmonicity in the system [29|. 



The expansion of the Compton profile in Eq. (|32| ) has a definite symmetry. 
The Hermite polynomial and the spherical harmonic each have parity (—1)', and 
hence J(— Q,— y) = J(Q,y). This is a symmetry which is obeyed by J(Q,y) 
if the impulse approximation is valid. Measured data, of course, will deviate 
from perfect symmetry. The procedure of fitting the symmetric expansion will 
ignore the antisymmetric components in the data, i.e. it will implicitly perform 
a symmetrization. This has, as a bonus, the effect of removing leading order 



corrections to the impulse approximation (see Sec. 2J3). For scattering from 
single crystals, there are likely to be additional symmetries in the directional 
Compton profile arising from the point group symmetry of the lattice site of the 
struck atom. In this case, the additional symmetry may be taken into account 
by expressing the expansion in terms of lattice harmonics [^] . 

As a final point, we discuss scaling of the data. A uniform scale factor yo is 
already included in Eq. ( |32|) . Choosing yo to be the width of the best Gaussian 
fit to the Compton profile matches the expansion to the data, and is likely to 
minimize the number of coefficients needed. However, this simple prescription 
is not adequate for strongly anisotropic data, where the width of the profile is 
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very different in different directions. Here we suggest a general procedure for 
anisotropic scaling. 



There does not appear to be a simple relationship between the inverse Radon 
transforms of two functions J(Q,y) differing by a linear coordinate transforma- 
tion, i.e. it does not appear useful to apply a linear transformation to the data. 
Instead, we ask the reverse question of how J(Q, y) changes when we trans- 
form the momentum distribution. Let us suppose, then, that the anisotropic 
distribution n(q) can be made reasonably isotropic by a linear scale change 

q-» q' = (qx/atx,qy/aty,q z /at z ), (35) 
for a suitable choice of orthogonal axes xyz. Defining n'(q) = n(q') and 

Q' = (a x Q x ,a y Q y ,a z Q z ), (36) 
we find the Radon transform of n' to be 

J'(Q, y) oc J n(q')%' • Q' - y) dq' = (l/Q')J(Q', y/Qf), (37) 

where Q' is the length of Q and Q' = Q'/Q', and an irrelevant constant factor 
has been omitted. The suggested procedure is then as follows. Given a mea- 
sured Compton profile, J(Q,y), let a x , a y and a z be its characteristic widths 
along the chosen axes. Compute the transformed profile J'(Q,y) according to 



Eq. (37), which amounts to a nonlinear coordinate transformation. After ap- 
plying the expansion technique described previously to reconstruct n'(q) from 
J', use 

ra(q) = n'(a x q x ,a y q y , a z q z ) (38) 

to find the required momentum distribution corresponding to the original J(Q, y). 
The transformed Compton profile J 1 is roughly isotropic. For example, if Q 
points along the x-axis, then J'(Q,y) oc J(Q,y/a x ), which has unit width in 
y, by construction. Thus J' is a function for which the Radon inversion step, 
using the expansion method, is efficient and stable. 



3.2 Anharmonic potentials 

The Compton profile for scattering from a harmonically bound atom is Gaus- 
sian, and a measurement of its width yields the mean kinetic energy. Essen- 
tially the same information might be obtainable by, for example, vibrational 
(infrared) spectroscopy. Where the neutron Compton scattering technique for 
bound atoms comes into its own, then, is in the potential to measure the an- 
harmonicity of the binding potential, in cases where it is legitimate to treat the 
struck atom in terms of single-particle motion in an effective potential. It is the 
only technique capable of measuring Born-Oppenheimer potentials directly. 

The extraction of the potential from the momentum distribution is simply a 
matter of inverting the Schrodinger equation: 

V(R) -E = (h 2 /2M)iJj'~ 1 'V 2 ip, (39) 



20 



-4 -2 2 4 

x/a 

Figure 3: Model wave functions for an atom in a double- well potential. Each is the 
sum of two Gaussian components of width a, centred a distance 2a apart. Solid line: 
a/a = 0.25; dashed line: a/a = 0.6; dotted line: a/a = 1. 



where ip(R) and ^(q) are Fourier transforms of each other, and the latter is 
obtained from n(q) = |-(/>(q)| 2 . The phase ambiguity in extracting the momen- 
tum space wave function from ra(q) is not a problem for inversion symmetric 
potentials, since ^>(q) can always be chosen to be real. 

It should be noted that the technique, at least in its simplest form, is re- 
stricted to systems in which all atoms of the species under study have identical 
chemical environments. If they do not, the Compton profile will be a superpo- 
sition of contributions from different environments. 

An extreme case of an anharmonic potential is one with two minima, such 
as might be expected for a proton in a hydrogen bond [28,29]. A useful model 



for illustrating the results expected for a double-well potential is to assume a 
ground state consisting of two shifted Gaussians of equal amplitude, 

= e -(*-a) 2 /^ 2 + e -(-+«) 2 /2a 2 ? (40 ) 

in one dimension |2S,2^,66|. This function is plotted in Fig. 3 for selected values 
of the ratio a/ a, and the corresponding potentials V(x) — E appear in Fig. 4. 
For the largest value oi a/a the potential consists of rather isolated harmonic 
wells, which gradually merge as a/ a is reduced. For a < a, the potential is 
more accurately described as a single harmonic well with a shallow "bump" in 
the centre, and the wave function is a single non-Gaussian peak. 

The Compton profile, J(y) = e~ a2y2 cos 2 ay, is plotted in Fig. 5 for an in- 
termediate value a/a = 0.6. It includes an oscillatory factor, representing 
interference between wave functions localized in the two wells, and the overall 
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Figure 4: Model potentials V(x) — E for each of the model wave functions in Fig. 3. 



shape is far from Gaussian. The number of oscillations in each half-width of the 
Gaussian envelope is roughly a/ a, so for weak anisotropy (a/a < 1) the devia- 
tions from Gaussian form are less pronounced. The coefficients in the Hermite 
expansion, Eq. (]32]), corresponding to this J(y) are A n oc (— a 2 /4a 2 ) n /(2n)\, 
which fall off rapidly even for intermediate a/ a. In the limit of large a/ a, the 
oscillations are so rapid that one observes only the average, which is Gaussian 
as one would expect for isolated harmonic well. 

This simple model may be extended to the case where the double well po- 
tential, as in a typical hydrogen bond, is not symmetric. It is found that the 
asymmetry supresses somewhat the oscillatory component in the wings of the 
Compton profile, which no longer goes to zero at its minima. 



4 Discussion 



The Compton profile in Fig. 5, corresponding to a model of an atom in a double- 
well potential, is far from Gaussian, which strongly suggests that anharmonicity 
of this kind should be experimentally measurable. Reiter p8| , p9[ has carried 
out an analysis of the practical limitations on measurements of anharmonicity, 
including the effect of sampling noise (but not of nonzero instrumental resolu- 
tion), concluding that extracting meaningful atomic potential energy functions, 
using the series expansion method to analyze data, is feasible. On the other 
hand, Sivia and Silver |B7|, in an analysis based on Bayesian probability theory, 
have shown that the reconstruction of the momentum distribution from the 
Compton profile is an intrinsically ill-posed problem, at least for studies of the 
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Figure 5: Compton profile corresponding to the two-Gaussian model wave function, 
with a/a = 0.6. 

condensate fraction in superfluid helium. A simple statement of the essence of 
their result is that widely different momentum distributions can yield the same 
Compton profile within experimental error. This does not, of course, imply 
that momentum distributions cannot be extracted. There are grounds to be- 
lieve that the inverse problem for strongly anisotropic data for single crystals 
is better behaved than for liquid 4 He: if the atomic motion is approximately 
decoupled into independent motion along three different axes, the resulting one- 
dimensional inversions are more stable than the three-dimensional inversion for 
isotropic data. Clearly, further study of the reconstruction problem for neutron 
Compton scattering from solids is desirable, to establish, for example, which 
features of the Compton profile are most sensitive to anharmonicity, the limits 
to accuracy of reconstruction, and optimum algorithms. 

Another question that deserves further investigation is the inclusion of cor- 
rections to the impulse approximation in the reconstruction procedure. This 
is a well-studied problem for scattering from superfluid helium [11], while for 



solids it is argued, as in Sec. 3A, that the corrections are negligible. Neverthe- 
less, we note that in Fig. 5 it would appear to be the tails of the profile which 
are most sensitive to the presence of anharmonicity in the potential, and even 
if corrections to the impulse approximation are small relative to the recoil peak 
amplitude, they may be significant in the tails. 

The past achievements of the neutron Compton scattering technique are im- 
pressive, ranging from the mature body of work on superfluid helium, to more 
recent applications in studies of quantum and classical liquids, and kinetic en- 
ergies, anisotropy and quantum effects in solids. Here we have summarized the 
theoretical background, and the prospects for a new development of the tech- 
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nique, the reconstruction of momentum distributions by direct inversion from 
the Compton profile. Taking into account theoretical work to date, there are 
grounds for being cautiously optimistic. We await the first direct experimental 
measurement of a three-dimensional atomic potential energy function in a solid. 
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